Modification of scattering lengths via magnetic dipole-dipole interactions 
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We propose a new mechanism for tuning an atomic s-wave scattering length. The effect is caused 
by virtual transitions between different Zeeman sublevels via magnetic dipole-dipole interactions. 
These transitions give rise to an effective potential, which, in contrast to standard magnetic interac- 
tions, has an isotropic component and thus affects s-wave collisions. Our numerical analysis shows 
that for ^"Cr the scattering length can be modified up to f 5 %. 
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Ultracold atomic collisions are essentially s-wave col- 
lisions. The s-wave scattering length a is the parame- 
ter that fully characterizes short-range interactions in an 
atomic Bosc-Einsten condensate (BEG) 0. Since it is 
deisrable to control the interactions in a BEG, tuning 
the scattering length has become an important issue. In 
bulk systems, a can be tuned by means of Feshbach reso- 
nance 5| . Additionally, low-dimensional systems exhibit 
confinement-induced resonance j^i^J. 

These methods do not affect the dipole-dipole long- 
range interactions. It is widely believed that short- 
range and dipole-dipole interactions act independently 
under usual experimental conditions. The influence of 
static magnetic dipole-dipole interactions on the proper- 
ties of a trapped BEG has been studied theoretically 
and experimentally . Laser- induced, finite- wavelength, 
dipole-dipole interactions were theoretically studied as 
well 0. The angular dependence of the static dipole- 
dipole potential is given by the second-order Legendre 
polynomial and thus does not directly affect s-wave col- 
lisions. Its indirect (mediated by the d-wave channel) 
influence changes significantly the s-wave atomic scat- 
tering length only if the interacting static electric dipole 
moments are induced by an external electric field of the 
order of 100 kV/cm which is a rather high value 
for an ultracold- atom experiment at present. The same 
mechanism is expected to be more efficient for highly- 
polarizable molecules . 

The conclusion that s-wave scattering is hardly af- 
fected by dipole-dipole interactions has been reached un- 
der the assumption that the dipole moments of the in- 
teracting atoms are given by the mean values of the cor- 
responding operators, while transitions to other Zeeman 
sublevels have been neglected. In the present Letter we 
demonstrate that these (virtual) transitions may lead to 
an appreciable change in the value of a. In this context, 
we discuss a remarkable general effect that has been over- 
looked thus far: the nonlocal collisional dynamics of ul- 
tracold atoms coupled by interaction between fluctuating 
magnetic dipoles. 

Since the magnetic moment of an atom in a weak mag- 
netic field is given by gF, g being the gyromagnetic factor 



and F the total angular momentum, the dipole-dipole in- 
teraction operator is 



VDD{r) = 



Hog^ FiF2-3(Fie,)(F2e, 



(1) 



where r is the interatomic separation vector, — 
r/r. Suppose the degeneracy of states with different z- 
projection of the momentum, = M, is lifted by an 
external magnetic field. Let the atoms be in the low- 
est energy state \F, M — —F). The excitation energy for 
the Zeeman sublevel |F, M = —F + 1) is denoted by /S.E. 
The operator couples the two-atom spin state 



= |F,M = -F)i|F,Af--F)2 



(2) 



to 



\(3) = {1/V2)(\F, M = ~F)i\F,M = -F + 1)2 

+ \F,M = -F+l)i\F,M = -F)2). (3) 

We adopt the two-channel approximation by represent- 
ing the wave functions of the colliding pair as = 
i/ja{j^)\a) + ■0/3(1") 1/3). Then the Schrodinger equation 
reads as 
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■V'a(r) 
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■V'/3(l-) 



{a\VDD{rmM^ 



(4) 



\VDD{rmM^) 



+ (/3|VDD(r)|«)Va(r), 



(5) 



Here m is the atomic mass, fi^k"^ /m is the collision energy, 
and ft^q^ — AE — h^k'^/m. Since we consider ultracold 
collisions, we have 



n^q^/m « AE. 



(6) 



The next step is to assume that AE is large enough 
to allow us to neglect the term (/3|yD£)(r)|/3)V'/3(r) m 
Eq. inj. Then the virtually excited state |/3) can be elim- 
inated using the negative-energy Green function of the 
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free motion G{-AE,r) = -e-«l'"l/(47r|r|) [T^: 
^^(r) = f d\'G{-AE,r-r'){P\VDD{r')\a)^Po.{r')■ 



Substitut ing Eq. into Eq. we obtain a non-local 
equation PJ for ipa- 

Va(i-) = V"Va(r) + (a|Vz)z?(r)|a)7/.„(r) 

m m 

+ (a|Vi)i,(r)|/3) j d'v'Gi-AE.v-v') 

x{l3\VDD{v')\a)i;^{v'). (8) 

On the assumption (verifiable a posteriori) that both G 
and (/9|VD£)|a) vary faster than ij^a with r', we can pull 
the wave function out of the integral and calculate the 
necessary matrix elements on expanding G'(— Ai?, r — r') 



as a series in spherical harmonics [l^- The resulting 
Schrodinger equation for a pair of colliding atoms with 
maximal projection of the angular momentum still re- 
tains the signature of nonlocal dynamics, by virtue of the 
spatial dependence of the effective potential on the con- 
volution of the Green function G(— Ai?, r — r') and the 
matrix element (/9| VD£i(r')|a). It has the form 



Va(r) = V^Va(r) 

m m 



2tp2 



4:TTr^ 

Qylmqg^F'^ 2 n- 2 



(1 - 3cos2 0)V'a(r) (9) 
cos^6sm^6W{qr)^Jc,{Y), 



where 9 is the angle between r and the quantization axis 
and 



W{qr) 



3 



sinh ^ 



2 cosh^ 



+ 



exp(-o(3r' + 3r' + r') 



exp{-qr)[3{qr) ^ 
3 



1 
qr 



sinh qr 



i{qr)-^ 
3 



q2j,2 



cosh qr 



(10) 



The function W{qr) is plotted in Fig. 1. It has two 
asymptotes: 



W{qr) « (qr)' 



qr > 1, 



(11) 



which corresponds to adiabatic elimination of in 
the limit wherein the approximation G{—AE, r — r') « 
—q'^S{r — r') holds (dashed line in Fig. 1), and 



W{qr) fa {6qry 



qr <t:l 



(12) 



(dot-dashed line in Fig. 1). 

To analyze the influence of the additional potential 
mediated by static magnetic dipole-dipole interactions, 
we average Eq. (|10(l over all angles and obtain in the 
ultracold-coUision limit (k 0) the following equation: 



J.2 ^j. 



dr 



^qll^W(qr) 
8r3 



where the quantity 

l„, = /iom.g2FV(4^?j') 



(13) 



(14) 



is the characteristic magnetic length. Equation H13() must 
be supplemented by the hard-sphere boundary condition 



V'Q(a) = 



(15) 




FIG. 1: Dimensionless function W vs. dimensionless distance 
qr. The two power-law asymptotes are shown (see details in 
the text). 



that accounts for the short-range (van der Waals) in- 
teractions. Solution of Eqs. H13() . I|15() yields the modi- 
fied scattering length a, which can be inferred from the 
r — > oo asymptotics of the solution that is propoportional 
to 1 — a/r. 

It is noteworthy that for g ^ the magnetically- 
induced isotropic potential varies as r^^ and does not 
depend on q. In this case Eqs. (|13|l . H15|) admit an ana- 
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FIG. 2; Modified scattering length for ^ Cr vs. frequency 
splitting A// = AE/{2nh). The two asymptotes of a in the 
small- and large-g limits are shown. 



lytical solution, 



4a 



so that the modified scattering length becomes 



V3?„,/(4a) 
tan[\/3Zm/(4a) 



(16) 



(17) 



In Fig. 2 we plot the results of numerical integration 
of Eq. (|13|l for ^°Cr. This isotope is chosen because 
chromium atoms have quite a large magnetic moment, 
gF = 6/iB, fiB being the Bohr magneton. On the other 
hand, the bare (unaffected by the effect under consider- 
ation) scattering length a for this isotope in the state of 
the maximum spin projection is of about 29 Bohr's radii, 
as can be concluded by mass-scaling arguments from the 
most recent data on ^^Cr scattering Il3|. Note, that the 
recent, more accurate measurements |l4| yield smaller 
values of the scattering length for ^^Cr and, con seq uently, 
for ^"Cr, than the previously reported data [l^. The 
value a = 1.54 nm is smaller than Im = 2.34 nm, and the 
difference between a and d becomes observable (up to 15 
%). By contrast, the scattering length for ^^Cr is three 
times larger than for ^°Cr. As a result, the influence of 
static magnetic dipole interactions in the case of ^^Cr is 
negligible. 

The discussed mechanism of the scattering length mod- 
ification will manifest itself via excitation of BEG oscil- 
lations by periodic modulation of the external uniform 
magnetic field. Such a method to reveal the scattering 
length modulation has been proposed with regard to Fes- 
hbach resonance |^^. Consider a BEG of '"^"Gr in the 
Thomas-Fermi regime [l^ trapped in an optical dipole 
harmonic potential, which is assumed, for the sake of 
simplicity, to be isotropic, loo being its fundamental fre- 
quency. The uniform, z-directed magnetic field oscillates 
in time as 



FIG. 3: Relative (dimensionless) modulation of the Thomas- 
Fermi BEC radius vs. dimensionless time under the action of 
the variable magnetic field Eq. (|THJ. (a) Bo = 0.08 mT, (b) 
Bo = 0.6 mT; Bi = 0.04 mT in both the cases. 



Its variation results in variation of the frequency splitting 
Ai' between the states \a) and The Zeeman shift in 
a weak field is given by Av/B = 87.7 MHz/mT. Thus 
we can infer the time variation of the scattering length 
using the dependence a(Aj/) plotted in Fig. 2. 

Let the magnetic field oscillations be switched on at f = 
0. The lowest monopole (angle-independent) mode of the 
BEC oscillations is described by the nonlinear equation 



(19) 



B{t) = Bo + Bi smujt. 



(18) 



where b = R/Rq is the ratio of the Thomas-Fermi ra- 
dius R of the oscillating BEG to the BEG radius at 
rest, R q = R\t= o = IhoiNdo/lhoV^^ |I3i = a\t=o, 
Iho = ^/h/JmuJo), N is the number of atoms in the BEC. 
For the numerical estimations we take u!o/{2tt) ^ 100 Hz, 
iV 10^, thus obtaining i?o ~ 6 fiui and the peak density 
of the BEC of about 3 • 10^'* cm"'^. The initial conditions 
to Eq. (0 are b\t=o = 1, &|t=o = 0. 

Choosing the frequency of the magnetic field modula- 
tion oj = V5wo, we ensure the resonance of the paramet- 
ric excitation with the lowest monopole mode in the lin- 
ear regime. In Fig. 3 we present the numerical solution of 
Eq. H19|) for two different choices of the magnetic field pa- 
rameters. While the amplitude Bi of the magnetic field 
oscillations is the same, the mean value, Bq, takes two 
distinct values. In Fig. 3(a) Bq corresponds to the steep- 
est derivative da/dB, so that the BEC excitations are 
excited most efficiently. For our numerical estimations 
we take u!o/{2Tr) ^ 100 Hz, N ~ 10^, thereby obtaining 
the Thomas-Fermi BEC radius Rq ^ 6 fim and the peak 
density of the BEC of about 3 • lO^"* cm^^. Thus for the 
parameters of Fig. 3(a), after 80 ms the amplitude of os- 
cillations of the Thomas-Fermi BEC radius approaches 
0.84 /im. This amplitude is approximately twice larger 
than the resonant wavelength of the optical ^6*3 <-^'' P3 
transition in chromium, thus making the BEC oscilla- 
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tions detectable by optical imaging. By contrast, the 
value of Bq in Fig. 3(b) corresponds to a relatively small 
da/dB and, therefore, to much less efficient excitation of 
the BEC. 

To conclude, we have proposed a new mechanism of 
s-wave scattering length tuning and modulation, differ- 
ent from both Feshbach-resonance and confinement- 
induced resonance methods 0| • It involves a hitherto 
unexplored aspect of ultracold-atom collisions, namely, 
their nonlocality mediated by fluctuations of the mag- 
netic dipole moments. Magnetic dipole-dipole interac- 
tion causes virtual transitions between different Zeeman 
sublevels of the atomic ground state. Such a coupling to 
the virtual state induces, additionally to van der Waals 
forces, an attractive isotropic interatomic potential. If 
the energy gap AE between the occupied and virtual 
spin states grows to infinity, we recover the "bare" (un- 
modified) scattering length value a, which is calculated 
from the van der Waals (non-magnetic) potential only. 
As AE vanishes, the renormalized scattering length ap- 
proaches the limit ()17|l that is fully determined by the 
ratio of the magnetic length 1^., Eq. (|14|l . to a. An ap- 
preciable effect is expected for Im/a ^ 1, which is the case 
for ^"Cr atoms, where we expect lowering of the modified 
scattering length by up to 15 % with respect to a. 
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